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ABSTRACT 
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^ I We use the number counts of X-ray selected normal galaxies to explore their 

evolution by combining the most recent wide-angle shallow and pencil-beam deep 

samples available. The differential X-ray number counts, dN/dS, for early and 

late-type normal galaxies are constructed separately and then compared with 

the predictions of the local X-ray luminosity function under different evolution 

>• I scenarios. The dN/dS of early type galaxies is consistent with no evolution out 

^r^\ to 2 ~ 0.5. For late-type galaxies our analysis suggests that it is the sources 

^ I with X-ray-to-optical flux ratio log fx/fopt > —2 that are evolving the fastest. 

O i Including these systems in the late-type galaxy sample yields evolution of the form 

)S^ I ~ (1 + z)'^^'^ out to z ~ 0.4. On the contrary late-type sources with log fx/ f opt < 



-2 are consistent with no evolution. This suggests that the log fx/ f opt > 



Qh' population comprises the most powerful and fast evolving starbursts at moderate 

O ' and high- 2;. We argue that although residual low-luminosity AGN contamination 

^ • may bias our results toward stronger evolution, this is unlikely to modify our 

^ ■ main conclusions. 
>'■ 

/\ ' Subject headings: Surveys - X-rays: galaxies - X-rays: general 

1. Introduction 

In the last few years surveys carried out by the XMM-Newton and Chandra missions 
provided the first X-ray selected samples of normal galaxies out to cosmologically interesting 
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redshifts (Hornschemeier et al. 2003; Bauer et al. 2004; Georgakakis et al. 2004; Norman 
et al. 2004; Georgantopoulos et al. 2005). This development has opened the opportunity 
to explore the evolution at X-ray wavelengths of the dominant population of the Universe: 
normal galaxies powered by stellar processes rather than accretion on a supermassive black 
hole. The unique feature of X-rays is that they are the only tool available to directly probe 
the X-ray binaries and the hot gas in galaxies, providing complementary information on their 
evolution compared to other wavelengths (e.g. Ghosh & White 2001). 

The evidence above has motivated a number of studies which use different approaches 
to constrain the evolution of X-ray selected normal galaxies. Stacking of optically selected 
spirals in the Chandra Deep Field North provided the first evidence for evolution of these 
systems at X-ray wavelengths (Brandt et al. 2001; Hornchemeier et al. 2002). A more direct 
method, which however requires the detection of individual systems, is to estimate the X-ray 
luminosity function at different redshift bins (e.g. Norman et al. 2004). A simpler but 
equally instructive technique uses the observed number counts in comparison with model 
predictions assuming different forms of evolution (e.g. Ranalli, Comastri & Setti 2005). 
The general consensus from the complementary methods above is that the X-ray evolution 
of normal galaxies to z ^ 1 — 1.5 is similar to that obtained from other wavelengths (e.g 
Hopkins 2004). 

In addition to the studies above, which mainly concentrate on deep surveys (e.g. Chan- 
dra Deep Fields) aiming to select high-z normal galaxies, there have also been efforts to 
identify such systems at low redshifts (Georgakakis et al. 2004; Georgantopoulos et al. 2005; 
Georgakakis et al. 2006). These samples are complementary to deep surveys, providing a 
tight anchor point at low- 2; that is essential for any conclusions on evolution. Despite the 
significant progress in the field however, there is still limited number of studies that exploit 
all the complementary low and high-z samples available to constrain the evolution of normal 
galaxies over as a wide redshift range as possible. 

In this paper we expand on previous results on the X-ray evolution of normal galaxies by 
combining the most up to date wide-angle shallow surveys (probing on average lower-z) with 
the latest compilation of normal galaxies in the Chandra Deep Field North and South (Bauer 
et al. 2004). These samples are used to construct the differential X-ray number counts of 
normal galaxies over about 4 dex in flux and to compare against different evolution models. 
We further improve on previous studies by splitting the sample into early and late-type 
galaxies to explore their evolution separately. This is particularly important since different 
galaxy types are expected to follow different evolution patterns. Throughout this paper we 
adopt Ho = 70kms^^ Mpc"^, Hm = 0.3 and Q\ = 0.7. 
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2. Sample selection 

In this paper we combine normal galaxy samples compiled from independent wide-area 
shallow and pencil-beam deep surveys spanning a wide range of fluxes, /x(0.5 — 2keV) ~ 
10-^3 _;L0-^^ erg s-^cm-2. 

At the bright fiux end we use galaxies identified in (i) the Needles in the Haystack 
Survey (NHS; Georgakakis et al. 2004; Georgantopoulos et al. 2005), (ii) the 1st release of 
the XMM-Newton Serendipitous Source Catalogue (IXMM) prodcued by the XMM-Newton 
Survey Science Centre (see Georgakakis et al. 2006 for the galaxy selection) and (iii) public 
Chandra observations, analysed as part of the XAssist project (Ptak & Griffiths 2003), that 
overlap with the Sloan Digital Sky Survey (SDSS; Schneider et. al. 2005). For the latter 
sample, full details on the field selection and normal galaxy identification will be presented in 
a future publication. The galaxy catalog for this sample is electronically available from the 
web (http://www.astro.noa.gr/xray/XAssist-SDSS.txt). In brief, sources are detected in the 
0.3-8 keV band of ACIS pointings with exposure times > 20 ks, to avoid observations that are 
too shallow to be useful for our purposes. Fluxes in the 0.5-2 keV band are estimated from 
the 0.3-8 keV count-rates in the XAssist catalogue assuming a power-law spectral energy 
distribution with F = 1.8. The sensitivity curve in the 0.5-2 keV band for these observations 
is constructed using methods described in Georgakakis et al. (2006). 

The total area of the combined bright sample (IXMM+NHS+XAssist) is about 15deg^ 
with a sensitivity curve shown in Figure 1. A total of 68 normal galaxy candidates are 
selected by searching for sources with log fx/fopt < —2, soft X-ray spectral properties to 
avoid type-II Seyferts and optical spectra (available for 55 of the 68 sources) consistent with 
stellar emission rather than AGN activity. For the log fx/ f opt calculation we adopt the 
definition of Stocke et al. (1991) assuming F = 1.8 and B — V = 0.8. 

At faint fiuxes, /x(0.5 - 2.0 keV) ^ 10"^^ - lO^^^ergs"^ cm'^ we use the CDF-North 
and South (Alexander et al. 2003; Giacconi et al. 2002) to identify normal galaxy candi- 
dates. We focus in particular, on a total of 104 sources (28 in CDF-S and 76 in CDF-N) 
classified normal galaxies by Bauer et al. (2004) on the basis of their X-ray and optical 
properties (e.g. hardness ratios. X-ray luminosity, optical emission lines), which appear con- 
sistent with stellar processes. This sample has a less stringent X-ray-to-optical fiux ratio 
cut, log fx/ f opt < —1, compared to the galaxies identified in the shallow surveys above. 
Because of evolution and/or differential X-ray /optical k-corrections, a more stringent X- 
ray-to-optical flux ratio cut (e.g. log fx/f opt = —2) is likely to affect more seriously these 
deeper samples which probe on average higher redshifts, resulting in incompleteness. On 
the contrary, the log fx/ f opt cut is likely to has little impact in the brighter lower-z samples 
described above (e.g. Georgantopoulos et al. 2005; Tzanavaris et al. 2006; see section 4). 
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3. Galaxy classification 

The normal galaxies identified in the surveys above are classified into early and late 
types by fitting galaxy templates to their broad band colours. For the X-ray bright samples 
presented here we use primarily the SDSS-DR4 ugriz photometry for spectral fitting. A 
total of 12 systems from the IXMM catalogue do not have SDSS observations available, in 
which case we rely on visual inspection of their optical spectra for classification. For the 
CDF-South we use the COMBO-17 photometry (e.g. Wolf et al. 2003) while for the CDF- 
North we also employ the SDSS-DR4 ugriz photometry. The exception are 7 CDF-North 
systems fainter than the SDSS limit, for which we use the multiwaveband observations of 
Barger et al. (2003) UBVRIzHK. 

The template SEDs are constructed by smoothly interpolating between the 4 galaxy 
types of Coleman, Wu & Widman (1980; E/SO, Sbc, Scd, Im) extended in the UV and IR 
wavelengths using the GISSEL98 code (e.g. Bruzual & Chariot 2003). This method provides 
a total of 60 SEDs from ellipticals (spectral classification 0) to extreme starbursts (spectral 
classification 60). Sources with best-fit SEDs in the range 0-25 are classified early type. This 
spectral range is in good agreement with the visual classification of the optical spectra (e.g. 
absorption lines), when available. 

Table 1 summarizes the galaxy types found in the four different surveys. It is evident 
that the bright IXMM, NHS and XAssist samples which include only low redshift galaxies 
{z < 0.2), consist of roughly equal numbers of early and late type galaxies. In contrast, 
the fainter CDF surveys consist primarily of late-type galaxies. This suggests that late and 
early type galaxies have different evolution with cosmic time. A quantitative analysis and a 
detailed discussion of this point are presented in the next section. 

Table 1: Number of normal galaxies in the samples used in this paper. 



sample 


early 


late 


all 


CDF (full) 
CDF (restricted) 

XAssist 
IXMM+NHS* 

total (full) 
total (restricted) 


35 
12 
13 
23 

71 

48 


69 
22 
9 
21 

99 

52 


104 
34 
22 
46 

172 
102 



*2 systems in the IXMM sample remain unclassified. See Georgakakis et al. (2006) for details. 
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Fig. 1. — Solid angle as a function of limiting flux in the 0.5-2 keV band for the X-ray surveys 
used in this paper to identify normal galaxies. The dotted (red) line is for the IXMM+NHS 
sample described by Georgakakis et al. (2006), the dashed (green) line corresponds to the 
Chandra serendipitous survey presented in this paper while, the continuous (blue) curve is 
for the combined CDF-North and South. The total sensitivity curve of these independent 
surveys is plotted with the (black) long-dashed line. 
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Fig. 2. — Differential normal galaxy counts in the 0.5-2 keV spectral band. The crosses are 
the combined sample of normal galaxy candidates including all systems from the CDF- North 
and South without applying any cutoff in log fx/ f opt- The continuous line is the prediction 
of the Georgakakis et al. (2006) XLF for both early and late type galaxies assuming the 
maximum likelihood evolution of the form oc (1 + z)^'^^. The open circles represent the 
subsample that only includes CDF galaxies with log fx/f opt < —2. The corresponding 
dashed line is the no evolution prediction {p = 0) of the Georgakakis et al. (2006) XLF. 
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4. The galaxy X-ray number counts 

In this section we discuss the number density of X-ray selected normal galaxies over 
4 decades of flux in the context of evolution models. We choose to use the differential 
galaxy counts, dN/dS, instead of the cumulative log A^ — log S", frequently employed in X-ray 
astronomy because of the independence of individual flux bins simplifying the estimation of 
errors and the interpretation of the results. 

Some contamination from low-luminosity AGN is expected in the ultra-deep Chandra 
fields used here, which comprise, on average, higher-z systems that are harder to study in 
detail. Therefore, the full Bauer et al. (2004) sample used here (total of 104 galaxies) 
likely represents close to the upper bound in the number density of X-ray selected normal 
galaxies. We also independently consider galaxies in the Bauer et al. (2004) sample with low 
X-ray-to-optical flux ratio, log fx/ f opt < —2. The AGN contamination in this restricted 
subsample is expected to be minimal, in the expense of possibly missing powerful starbursts 
or massive ellipticals, suggested to occupy the log fx/ f opt > —2 region of parameter space 
(e.g. Alexander et al. 2002; Bauer et al. 2002). More recently, Tzanavaris et al. (2006) 
cross-correlated the 2dF Galaxy Redshift Survey with public XMM-Newton and Chandra 
observations and found a non-negligible fraction normal galaxies at z ^ 0.1, mostly early 
types however, with log fx/f opt > —2. The evidence above suggests that the Bauer et 
al. (2004) CDF-North and South restricted subsample with log fx/ f opt < —2 is likely to 
represent the lower bound in the number density of normal galaxies. 

Figure 2 presents the 0.5-2 keV differential counts for all normal galaxy candidates in the 
flux range /x(0.5— 2.0 keV) ~ 10^^^ — 10^^^ erg s~^ cm~^ identifled in the surveys described in 
section 2. The dN/dS in this flgure is constructed using the sensitivity curve shown in Figure 
1. We consider separately the two samples that include (i) all the Bauer et al. (2004) galaxies 
and (ii) only those systems in the CDF-North and South with log fx/ f opt < —2. The former 
is hereafter referred to as the full sample while, the latter is dubbed the restricted sample. 
We flt the unbinned data for these two galaxy samples with a power-law using the parametric 
Maximum Likelihood method. This yields slopes of —2.47 ± 0.04 and —2.22 ± 0.05 for the 
full and restricted samples respectively. Using the full sample and adopting the combined 
early and late-type X-ray luminosity function (XLF) of Georgakakis et al. (2006) assuming 
luminosity evolution of the form (1 -|- zY we flnd p = 1.92 ± 0.35. This is consistent with 
the evolution derived by Norman et al. (2004) and Ranalli et al. (2005). The dN/dS of 
the restricted sample with log fx/ fopt < —2 is best flt with negative values of p, implying 
no-evolution. The evidence above suggests that in the case of the total sample (both early 
and late types) the evolution is driven by the log fx/ fopt > —2 systems. 

In addition to the total galaxy counts we also consider early and late-type galaxy sub- 
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Fig. 3. — Differential counts for early-type normal galaxies in the 0.5-2 keV spectral band. 
The crosses are evolved galaxy candidates, including all systems from the CDF-North and 
South without applying any cutoff in log fx/ fopt- The open circles represent the subsample 
that only includes CDF early-type galaxies with log fx/ fopt < —2. The continuous line is 
the prediction of the Georgakakis et al. (2006) XLF for early-type galaxies assuming no 
evolution. 
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Fig. 4. — Differential counts for late-type normal galaxies in the 0.5-2 keV spectral band. 
The symbols are as in Figure 3. The continuous and dashed lines are the predictions of the 
Georgakakis et al. (2006) late-type galaxy XLF for luminosity evolution oc (1 + z)^''^ (i.e. 
the maximum likelihood fit to the full subsample) and no-evolution respectively. 
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samples independently. Figure 3 plots the dN/dS for evolved galaxies according to the 
classification presented in section 3. The full and restricted subsamples for early- type galax- 
ies are plotted separately, delineating the upper and lower X-ray number density limit of 
evolved galaxies. In the case of the full subsample, using the parametric maximum likeli- 
hood method and assuming luminosity evolution of the form oc (1 + zY for the Georgakakis 
et al. (2006) early-type XLF we constrain p = O.SQlg'gg consistent with no evolution. The re- 
stricted subsample requires negative values of p implying either incompleteness of the CDFs 
due to the log fx/ f opt < ~2 cutoff or uncertainties in the early/late type classification. Nev- 
ertheless the evidence above suggests that optically classified early-type galaxies show little 
(if any) evolution at X-ray wavelengths, at least out to z ~ 0.5, the median redshift of 
evolved galaxies in the CDFs. 

In Figure 4 we plot the differential counts for late-type galaxies. These are compared 
with the dN/dS estimated by integrating the late-type galaxy XLF of Georgakakis et al. 
(2006) assuming different scenarios for the evolution of this population. The restricted 
subsample, comprising systems with log fx/fopt < —2, appears in agreement with the no 
evolution model. Indeed, a maximum likelihood fit to the data, parameterising the evolution 
as oc (1 + zY for the Georgakakis et al. (2006) late-type XLF, yields p = 0.4:1~^q^I consistent 
with no evolution. On the contrary, a similar analysis for the full subsample gives p = 2.70± 
0.35, suggesting evolution out to z ~ 0.4, the median redshift of the late-type subsample. 
This suggests that it is the log fx/fopt > —2 galaxy population that evolves the faster. 
Indeed, these sources are suggested to comprise a large fraction of starburst galaxies expected 
to show most dramatic evolution with redshift (e.g. Alexander et al. 2002). The evolution 
derived here for the late-type sample is comparable within the uncertainties to this of the 
full sample. This is not surprising as most of the galaxies at high redshift, which drive 
the evolution, are late-types. The derived evolution at X-ray wavelengths is consistent 
with that estimated for star-forming galaxy samples selected at other wavelengths (e.g. 
Hopkins 2004). Therefore, the late-type X-ray population studied here closely maps the 
star-formation history of the Universe, at least out to the median redshift of the sample, 
z ~ 0.4. The evolution rate derived here, p = 2.70 ± 0.35, is consistent with those models of 
Ghosh & White (2001) that adopt a slow evolution timescale for the low- mass X-ray binary 
population, > 1 Gyr. Shorter timescales predict exponents much higher than that derived 
here, almost independent of the star-formation history of the Universe at higher redshift. 

Finally, we caution that low-luminosity AGN contamination may be present in our 
sample. This problem is likely to be more severe for higher-^; systems, where weak AGN 
optical signatures are harder to identify against the stellar continuum (e.g. Severgnini et al. 
2003). Although, the presence of low-luminosity AGN may bias our results to higher values 
for the evolution parameter, p, we argue that this is unlikely to modify the main conclusions 
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presented here. The X-ray spectra of many low luminosity AGN are indeed, dominated by 
thermal emission below about 2keV (Levenson, Weaver & Heckman 2001; Terashima et al. 
2002), suggesting stellar processes. Even if the late-type galaxy sample used in this paper has 
some level of residual AGN contamination, this is likely to have only a minor contribution 
to the 0.5-2 keV, employed here. 
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